Superconductivity in the t-J model 
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A comparison of microscopic theories of superconductivity in the limit of strong electron correlations 
is presented. The results for the two-dimensional t-J model obtained within the projection technique 
for the Green functions in terms of the Hubbard operators and the slave-fermion representation for 
the RVB state are considered. It is argued that the latter approach resulting in the odd-symmetry 
p-wave pairing for fermions is inadequate. 



I. INTRODUCTION 

A mechanism of high-temperature superconductivity 
in cuprates is still unresolved since strong electron cor- 
relations in copper-oxygen planes prevent from applying 
well established methods of band structure calculations 
developed for conventional metals. An important role 
of electron correlations in cuprates was initially stressed 
by Anderson 1 who proposed to consider them within the 
framework of the Hubbard model or the so-called t-J 
model which follows from the Hubbard model in the limit 
of strong correlations. To study these models, a lot of 
numerical work have been done 2 ' 3 though the obtained 
results are still controversial. For instance, a robust d- 
wave pairing was observed for the t-J model 4 , while a 
long-range order was not found in the original Hubbard 
model. 5 

In analytical approache mostly a mean-field theory 
(MFA) was applied in studies of the Hubbard model or 
the t-J model. The resonating valence bond (RVB) state 
in the t-J model was proposed by Baskaran et al. 6 where 
superconductivity was obtained as a result of spin corre- 
lations induced by the superexchange interaction. Sim- 
ilar results were found by Cyrot 7 for superconducting 
pairing mediated by superexchange interaction. To over- 
come the problem of strong correlations in the t-J model 
and nonfermionic commutation relations for the physical 
electron operators a number of auxiliary field represen- 
tations was proposed (see, e.g., 8-32 ). However, in these 
methods usually a spin-charge separation is assumed for 
spinon and holon fields which violates rigorous restric- 
tions imposed by nonfermionic commutation relations as 
the "no double occupancy" constraint which may result 
in unphysical conclusions. 

A reliable analytical approach to deal with strong cor- 
relations in the Hubbard or the t-J model is based on 
the Green function methods in terms of the Hubbard 
operators which rigorously preserve the non-fermionic 
commutation relations 33 . Here we may mention our re- 
sults based on the Mori-type projection technique for the 
Green functions 34-36 and the diagram technique calcula- 
tions by Izyumov et al. 37,38 . These approaches enable one 



to go beyond MFA by taking into account self-energy cor- 
rections. For instance, a numerical solution of the Dyson 
equation in 36 revealed a non Fermi-liquid behavior in the 
normal state at low doping and the d-wave supercon- 
ductivity mediated by the exchange and spin-fluctuation 
pairing. In the recent paper 39 a microscopical theory of 
superconductivity in Cu02 layer within the effective two- 
band p-d Hubbard model in the strong correlation limit 
was developed. It has been proved that the MFA for the 
Hubbard model results in the antiferromagnetic exchange 
d-wave pairing which is equivalent to pairing observed in 
the t-J model in MFA. 

In the present paper I compare the results for the t-J 
model obtained within the Green function (GF) method 
in terms of the Hubbard operators 34-36 and by apply- 
ing the slave-fermion hard-boson representation 2 for the 
RVB state. It will be shown that the latter approach 
results in an odd-symmetry p-wave pairing for the spin- 
less fermions as in the path-integral representation in 32 
which contradicts to known numerical and analytical cal- 
culations. It casts doubts on the validity of spin-charge 
separation approach in studying superconducting pairing 
in the t-J model. 

In the next Section we present briefly the results of the 
projection technique for the GF 36 for the t-J model. In 
Sec. 3 superconducting pairing within the slave-fermion 
representation for the Hubbard operators is considered. 
Results and discussions are given in Sec. 4. Concluding 
remarks are in Sec. 5. 



II. GREEN FUNCTION METHOD 

A. Dyson Equation for the t-J Model 

In the present Section we consider superconducting 
pairing in the t-J model by applying the GF tech- 
nique 34-36 . The t-J model in the standard notation 1,42 
reads: 



H t -j = -* X) + 



(1) 



i 



where a\ a = a^(l — nj_ CT ) are projected operators for 
physical electrons, rij = ^ CT af a a i(J is the number opera- 
tor and Sf = (1/2)X) S)S ' ^ct" s ,(W are spin-1/2 opera- 
tors. Here t is effective transfer integral and J is antifcr- 
romagnetic exchange energy for a pair of nearest neighbor 
spins, (ij), i> j . 

To take into account on a rigorous basis the projected 
character of electron operators we employ the Hubbard 
operator technique. The Hubbard operators (HO) are 
defined as X" 13 = \i,a){i,(3\ for three possible states 
\i, a) at a lattice site i: for an empty site \i, 0) and for a 
singly occupied site \i,cr) by an electron with spin a/2 
(er = ±1, a = —a). They obey the completeness relation 



(2) 



which rigorously preserves the constraint of no double 
occupancy. The spin and density operators in Eq.(l) are 
expressed by HO as 



(3) 



The HO obey the following multiplication rule 
X^X^ 5 = Sp-yX" 6 and commutation relations 



i j j 



% faxf ± 5 5a xf) , (4) 



where the upper sign stands for the the Fermi-like HO 
(as, e. g., X® a ) and the lower sign for the Bose-like op- 
erators (as the spin and number operators in Eq. (3) ). 

By using the Hubbard operator representation we write 
the Hamiltonian of the t-J model (1) in a more general 
form: 



H, 



t-J 



<7<7 \ 

3 ) 



(5) 



The electron hopping energy for the nearest neighbors, 
tij = t, and the second neighbors, t^ = t , on a 2D 
square lattice, and the exchange interaction J,j = J for 
the nearest neighbors 43 can be considered as independent 
parameters if, starting from a more realistic for copper 
oxides three-band p-d model 44 , we reduce it to the t-J 
model 42 . In that case the parameters t, t' and J can be 
evaluated in terms of the original parameters of the p-d 
model (see, e.g., 45 ' 46 ). We introduced also the chemical 
potential /z which can be calculated from the equation 
for the average number of electrons 



n 



(6) 



To discuss the superconducting pairing within the 
model (5) we introduce the Nambu notation for HO: 



X 0a 

X 



and consider the matrix GF 



GiiAt t') = «<Mi)|*t (f)» = ( % fi° ) (7) 

where Zubarev's notation for the anticommutator GF is 
used 47 . 

By differentiating the GF (7) over the time t we get 
for the Fourier component the following equation 



u)G ija {w) = 5 lJ Q a + {{Z itT | * ja )] 

Qa 



where Z ia = \^ ia ,H\ , Q a 







(8) 



with Q a = 



+ Xf a ). Since we consider a spin-singlet state the 
correlation function Q a = Q = 1 — n/2 depends only on 
the average number of electrons (6). 

Now, we project the many-particle GF in (8) on the 
original single-electron GF 



+ ((zt rr) I 



jet Hi 



(9) 



where the irreducible iirr) part of the many-particle op- 
erator Zi a is defined by the equation 

{{Z?; r \*%}) = {{Z^-Y.Eu^i,, *+}> =0, 



which results in the definition of the frequency matrix 

E ij<7 = {{\* ia ,H],*+ a })Q-' . (10) 

Now we can introduce the zero-order GF in the general- 
ized MFA which is given by the frequency matrix (10) 



G° ija (Lj) = Q{LJT S ij -E ija }- 1 . 



(11) 



To derive the Dyson equation for the single-electron 
GF (7) we write down an equation of motion for the irre- 
ducible part of the GF in (9) with respect to the second 
time t' for the right-hand side operator ^~j a (t'). Then 
performing the same projection procedure as in Eq. (9) 
we obtain the Dyson equation for the GF in the form 

G ija (w) = G°- CT M +E^LM EfcfcrM G lja (w), (12) 
u 

where the self-energy operator Ejy^w) is defined by the 
equation 

f ij(T (u) = £ija(u) + E^( w ) Gkiai^) Ti ja (cj) . (13) 

kl 

Here the scattering matrix is given by 



2 



f ij ^) = Q- 1 ({zt ) \zS r) % Q- 1 - 



(14) 



From Eq. (13) it follows that the self-energy operator is 
given by the irreducible part of the scattering matrix (14) 
that has no parts connected by the single zero-order GF 
(11): 

t ija (w) = Q-'az^ i ^r )+ »r } q- 1 ■ (15) 

Equations (11), (12) and (15) give an exact representa- 
tion for the single-electron GF (7). However, to solve the 
self-consistent system of equations one has to introduce 
an approximation for the many-particle GF in the self- 
energy matrix (15) which describes inelastic scattering of 
electrons on spin and charge fluctuations. 



B. Self-Consistent Equations 

Here we derive a self-consistent system of equations in 
MFA. To calculate the frequency matrix (10) wc use the 
equation of motion for the HO: 

l.u' 



l,a' 



where we introduced the operator 

= (1 - ^ + 5?)^ + 5f W (17) 

The Bose-like operator (17) describes electron scattering 
on spin and charge fluctuations caused by the kinematic 
interaction (the first term in (16)) and by the exchange 
spin-spin interaction (the second term in (16)). 

By performing commutations in (10) we get for the 
normal and the anomalous parts of the frequency matrix: 

En = - M % + % J2i^( x r° x ? a )/Q 



+ ^Ju(Q-i + xu/Q)} 
UAQ + xVi/Q) \ja{xfx^)/Q (is) 



En = A ija = % ]T t u {XfX^ X°°Xf)IQ 



- -Ja (X^Xf - Xfxf)/Q . (19) 

Here in calculation of the correlation function for the 
normal component of the frequency matrix: 



((l-- n . + S ?)(l-- nj + S]) + S?SJ) 
= ((1 - \n )) + (SjSj) = Q 2 + X % , 

we introduce the charge- and spin-fluctuation correlation 
functions 

X% = \(5n i 6n j ) + (S i S J ) , 

with Srii = rii — (rii). Further we neglect the charge 
fluctuations, (SniSrij) ~ 0, but take into account spin 
correlations given by spin correlation functions for the 
nearest (xis) and the next-nearest (X2s) neighbor lattice 
sites: 



Xls ~ (SiSj +ai ) , X2s — (SiS i+a2 ), 



(20) 



where a\ — (±a x ,±a y ) and a 2 = ±(a x ± a y ). In a para- 
magnetic state assumed here they depend only on the 
distance between the lattice sites. 
In the k-representation for the GF 



wc get for the zero-order GF (11): 

G ( °Hk,uj) =Q{ujt - (e k - m)t 3 - A^}- 1 



= Q 



( £k -M) 2 -|A£ 



(21) 



where To, fi, 73 are the Pauli matrices. The quasipar- 
ticle energy e k and the renormalized chemical potential 
fi = fj, — 5fi in the MFA are defined by the frequency 
matrix (18) 



-t(k) 



2J 



^ 7 (k-q)iV ( 



^ = jf E - 2 W 2 - Xu/Q) 



(22) 



(23) 



where J(q) = 4J 7 (q), t(k) = 4tj(k) + 4i' 7 '(k) 
with 7 (k) = (l/2)(cosa x q x + cosa y q y ) , 7 '(k) = 
cos a x q x cos a y q y , while the renormalized hopping inte- 
gral is given by 



t(k) =4tj(k)Q(l + xis/Q 2 ) 
+ U>i{k)Q{l + X28 /Q 2 ). 



(24) 



The average number of electrons in Eqs. (22), (23) in the 
k-representation is written in the form: 

n K * = (X^°X^) = QN ka . 

It should be pointed out that the renormalization of the 
hopping parameter (24) caused by the spin correlation 
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functions (20) are essential at low doping when short- 
range antifcrromagnetic correlations are strong. For in- 
stance, for hole doping S = 1 — n ~ 0.05 and Q — 
(1 + 8)/2 ~ 0.53 the correlation functions are estimated 
as 36 : xis — —0.3, x?s — 0.2, which results in complete 
suppression of the nearest neighbors hopping, while the 
next-neighbor hopping is quite large: 

t ef/ =tQ[l + X is/Q 2 } ^0, 

t' eff =t' Q[l + X 2s/Q 2 ] ^0.9i'. 

For large doping the antifcrromagnetic correlations are 
suppressed and the nearest neighbor hopping prevails: 

teff /t'eff - t/t> » 1 . 

The superconducting gap A£ in Eq. (21) is defined 
by the anomalous component of the the frequency ma- 
trix (19): 

A £ = E - l J ( k 1] ( X %K a )> (25) 
q 

There are two contributions in Eq. (25) given by the k- 
independent kinematic interaction t(q) and the exchange 
interaction J(k — q). The kinematic interaction gives 
no contribution to the <i-wave pairing in MFA, Eq. (25) 
(see 34 ), and we disregard it in the following equations 48 . 
The anomalous correlation function in Eq. (25) can be 
easily calculated from the anomalous part of the GF (21): 

(A * X" a ) = -Q-^ 3 - tanh ^ , (26) 
\ - q q / 2T , v ) 

which results in the BCS-type equation for the gap func- 
tion: 

A£ = — J(k - q)-^2- tanh — , ( 2 7) 
k N ^ y 2E Q 2T ' v ' 

q 

where E k = [ (e k - /i) 2 + | A£ | 2 ]V 2 is the quasiparti- 
cle energy in the superconducting state. As was proved 
in 39 , the retardation effects for the exchange interaction 
are negligible and therefore there is no restriction in in- 
tegration over the energy in Eq. (27). It means that all 
electrons in the conduction band participate in the su- 
perconducting pairing contrary to the BCS equation for 
the electron-phonon model where the energy integration 
and pairing are restricted to a narrow energy shell of the 
order of the phonon energy close to the Fermi energy. 

The equation (27) is identical to the results in MFA of 
the diagram technique 37 , while the gap equation obtained 
in 40 has an additional factor Q = (1 — n/2) which is spu- 
rious. This factor appears if we apply a simple decoupling 
procedure in the equation of motion (16) for GF instead 
of the projection technique given by Eq. (10). Writing 
the bosonic operators in the exchange interaction as a 
product of two fermionic operators: Xf a = Xf°X^ a 
and performing decoupling of the fermionic operators: 



{{Xf°X^Xf - XfX^ s X^ | Xf)) u ~ 
(X°°Xf XfXf) «Xf° | Xf)) u , (28) 

we obtain the same equation (27) for the gap function but 
with the additional Q-factor at the right hand side. In 
the decoupling we miss the normalization factor for the 
correlation functions in the denominator of the frequency 
matrix, Eq. (10), which cancels out with Q in the numer- 
ator of the corresponding GF (21). So, a rigorous way to 
apply MFA with a proper account of projected charac- 
ter of HO is to use the projection technique as discussed 
above. 

The self-energy contribution (15) in the second order of 
the the kinematic interaction is considered in 36 while it is 
omitted in 37 . As discussed in 36 ' 39 , it mediates the spin- 
fluctuation pairing and results in finite life-time effects for 
the quasiparticle spectrum giving rise to an incoherent 
contribution to the single-particle density of states. Here 
we shall not discuss further these self-energy effects since 
for comparison the GF approach with the slave-fermion 
technique it will be sufficient to consider only MFA for 
the gap equation (27). 

III. SLAVE-FERMION APPROACH 
A. Slave-fermion representation 

A number of auxiliary field representations has been 
proposed so far (see, e.g., 8-32 ). In the slave-boson 
method 8-10 the projected electron operator is expressed 
as a product of the auxiliary Bose field for the charge 
degree of freedom (holon) and the Fermi field for the 
spin degree of freedom (spinon). The main problem in 
this approach is the so called constraint imposed by pro- 
jected character of electronic operators in the t-J model 
which prohibits double occupancy of any lattice site. To 
treat the constraint a site-dependent Lagrange multiplier 
is introduced. However, to solve the problem the MFA is 
usually applied and the Lagrange multiplier is taken to 
be independent on the lattice site so the local constraint 
is replaced by a global one with uncontrollable conse- 
quences. In the slave-fermion method the charge degree 
of freedom is represented by spinless fermion operators, 
while to describe the spin degree of freedom the Bose 
field (Schwingcr bosons 13-22 or spin operators 23-27 ) are 
used. The Schwinger boson representation though be- 
ing physically meaningful for the Heisenberg model 11 ' 12 
gives poor results for the doped case: the antifcrromag- 
netic ground state persists up to very high doping and 
it does not reproduce the large Fermi surface as in the 
slave-boson method. In the slave fermion and spin oper- 
ator representation the magnetic properties of the model 
are described in a more reliable way (see, e.g., 25 ' 29 ' 31 ). 

Below we consider the slave-fermion hard-core (CP 1 ) 
boson representation proposed in 27 and later on cm- 
ployed in 28 -31 to investigate different physical proper- 
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ties of cuprates within the t-J model. It has some ad- 
vantage since the constraint of no double occupancy can 
be fulfilled without introducing the Lagrange multiplier. 
To decouple the charge and spin degrees of freedom for 
physical electrons the HO in the theory 27 is represented 
as a product of a spinless fermion hf for the charge de- 
gree of freedom (holon) and a hard-core boson bi a for 
the spin degree of freedom (spinon): 

X^ = hfb itr , Xr° = h i bf, (29) 

which have the following commutation relations: 

hfhj + hjhf = 5 itj , b ia b+ a - b+ a b ia = S itj (l - 2b+b ia ). 

The hard-core bosons are Pauli operators which commute 
on different lattice sites and anticommute on the same 
lattice site prohibiting double occupancy. The Pauli op- 
erators can be also represented by the spin-lowering S[~ 
and spin-raising operators for spin- 1/2 : 

6+ = S+, 6, T = Sr, or b+ = Sf, b ta = Sf. 

The on-site electron local constraint 

Y x?xy = y xr - hh+ y ^ 

a a (7 

= \Kl = 1 - h+hi < 1 (30) 

is satisfied here since for the Pauli operators at any lattice 
site we have the equation 

Y h t bi« = + Sr S + = 1, (31) 

a 

and the spinless holon number n\ h ^ = h^h i can be equal 
to 1 or 0. 

However, the spin-charge separation imposed by the 
representation (29) results in extra degrees of freedom: 
a spin 1/2 is assigned to any lattice site including an 
empty site, while in the HO representation, Eq. (2), we 
have only 3 states: an empty state and a filled state 
with spin ±1/2. To cure this defect one should introduce 
a projection operator to exclude the unphysical states 27 
(The same applies to other slave-fermion spin operator 
representations as in 23 ~ 26 ). Otherwise the commutation 
relations for the original HO, Eq. (4), and their represen- 
tation, Eq. (29), will give different results. For instance, 
for the physical electrons which are described by HO we 
have 

\ A i i i S — A i + A i — 1 — A i ) 

and Y < W- X i°}) = 2 - K> = 1 + 5 > (32) 

(7 

where the hole doping concentration 5 = 1 — n . If we 
use the representation (29) then we can write the com- 
mutation relations as 



{X^,Xf} = {hfb i<T ,h i bi} 

= h+ hi + {\-2 h+ hi)b+ a b la 

and ^{/ l +6 iCT ,/i i 6+} = l, (33) 

where we have used Eq. (31). For the average value in 
Eq. (33) we get, respectively, 

Y({ h th*,h i b+}) = 1, (34) 

(7 

which contradicts to the rigorous result for HO, Eq. (32). 
For the average number of electrons in the representa- 
tion (29) by using the definition (6) we can write 

n = { ni )=Y(XfX° a ) 

a 

= I>i = (fHht) = 1-8, (35) 

which coincides with the definition (6) if we take for the 
hole number operator the following definition: X®° = 
Kl hi. However, this definition is not unique. For in- 
stance we can write: X™ = X^Xf = XfXf , which 
results in the equation: = (l/2)hfhi if we use the 
representation (29) and the conditional). So the dou- 
ble counting of empty sites results in controversial equa- 
tion for an average number of electrons which is therefore 
valid only with accuracy of ±5. 

B. Mean- field approximation 

Let us consider the resonating valence bond (RVB) 
state in the original Hamiltonian (5) as proposed by 
Baskaran et al. 6 . For this we should write the Bosc- 
like spin operators in the exchange energy in H t -j (5) 
as a product of two single-particle Fermi-like operators: 
Xf 7 — Xf®X® a and introduce the singlet operators 

bl^b^^X^Xl-X^X^). (36) 

Then using MFA for the singlet operators in the exchange 
interaction of the t-J model we get the RVB effective 
Hamiltonian: 

H ' = \ E MK-x; + xt+xj-) = -\Y -W>:, 

^4E J « {Bpij+btB.j-lB^ 2 ), (37) 

where we introduced the RVB order parameter: 
B<P ee Bij = (6 y ) = ^=(X9+X°- V t °-V°+). (38) 

Here we should point out that in MFA (37) a decoupling 
of the Hubbard operators on the same lattice site is used: 
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(X+°X7°) X° + X°~ , (39) 

which is not unique and results in uncontrollable approxi- 
mation since the Hubbard operators obey the multiplica- 
tion rule: X^X? 7 = A" 7 , and any intermediate state 
(3 can be used in the decoupling. 

To obtain a self-consistent equation for the order pa- 
rameter we assume a spin-charge separation as it is usu- 
ally done in the slave-particle methods by applying a de- 
coupling of spinon and holon degrees of freedom intro- 
duced in Eq. (29): 

Bij = -^={h+h ] h+b n - h+bHh+bn) 

~ ^ h P^ b ^ b n - Wit) = F tj (40) 
K = -j={htbi,hp n - h+bHhpn) * hfhpij. (41) 

Within these approximations we obtain the effective 
Hamiltonian for holons 

~ \ E J H « h i &*> h t K + h i Hhf (42) 

where the effective hopping parameter and the exchange 
energy are given by 

Uj = Uj (6+ b n + b+ bji) = Uj (S+ Sr + S r <?+) 

Jij = Jij \<Pij\ 2 = Jij \ I (bit b Jl - b n bjt)\ 2 . (43) 

To obtain a phase diagram for the RVB order parame- 
ter = F^j ipij as a function of temperature T and 
hole doping 8 one should solve a self-consistent system of 
equations for the both order parameters, holon Fij and 
spinon ifij ones. Here we consider only equations for 
the holon order parameter by suggesting that there ex- 
ists a region of (T, 8) where the spinon order parameter 
is nonzero. 

Introducing k-vector representation for the correlation 
functions: 

i^H« = ^£eW))F(k) 

k 

= ^£sin(k(i-j))<fc k ft_ k >, (44) 

k 

<PH = ^= (6iT bji b tl b n ) = ^]T jW-M p(k) 
= |^sin(k(i-j))-L(6 kT 6_ ki ), (45) 
the Hamiltonian (42) can be written in k-space as 



k 

-i£{A(k)^ k / i + + A+(k)/ lk h- k ], (46) 
1 k 

where the chemical potential for holons ^ = — fi and 
the holon spectrum e(k) according to Eq. (43) is written 
in the form (cf. Eq. (24)): 

e(k)=t(k)=4t 7 (k)2<S+Sr Hoi ) 

+ 4t' 7 '(k)2(S+S- Q2 >, (47) 

Here we assumed that the spin correlation functions have 
s-wave symmetry, e.L, Xi,i±a x = Xi,i±a y • We also intro- 
duced the holon gap function defined by the equation 

A(k) = - A(-k) = 1 £ J(k - q) (h* fc_ q ) . (48) 
q 

In Eqs. (47), (48) we used notation e(k), A(k) for the 
holon spectrum and the holon gap to distinguish them 
from that ones of physical electrons in Sec. 2.2. 

The normal and anomalous correlation functions can 
be easily calculated for the BCS-type Hamiltonian (46): 

w^-K 1 - £j ^r t - h ^)- (49 ' 

F(q) = (ft q />_,) = ^ tanh (50) 

where the quasipartical spectrum E(k) = [(e(k) — 
/i/j) 2 + | A(k) | 2 ] 1/2 . From these equations follow the self- 
consistent equations for the gap function and the average 
number of holes 8 which defines the holon chemical po- 
tential ^ 

q 

q 

= lylfl- £(q) -^Vanh^M. (52) 

These equations are identical to the results obtained in 32 
where the path-integral representation for the t-J model 
was applied and an effective BCS-type Hamiltonian anal- 
ogous to Eq. (46) was derived. Starting from the MFA 
for RVB order parameter as in Eq. (37), the authors in- 
troduced MFA for the spinon auxiliary field that resulted 
in the spin-charge separation as in our Eq. (41). There- 
fore, we have proved that the results of the path-integral 
representation for the t-J employed in Ref. 32 are equiva- 
lent to the MFA for the slave fermion - hard-core boson 
approach considered in this section. 
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C. Holon Green functions 



IV. RESULTS AND DISCUSSION 



To avoid uncontrollable approximation caused by the 
decoupling of the Hubbard operators on the same lattice 
site in Eq. (39) used in MFA for RVB state in the Hamil- 
tonian (37), in the present section we consider projection 
technique for the holon GF for the spinon-holon model. 
By using the spinlcss fermion hard-core boson represen- 
tation (29) we write the Hamiltonian of the t-J model (5) 
as follows: 



H, 



t-J 



- ji hj K\ — ^ Jij hi hj hj h\ 

x (5+ S; Sr S+)(S+ 5- - Sr <?+). (53) 

To obtain equation for the holon GF we apply the pro- 
jection technique described in Sec. 2.1. By introducing 
the matrix GF for holon operators: 

rh (f _ t ,,( «M*)l>tf(t')» «M*)IM*')» 
ij[ ' \{{htm+(t'))) «/>+(*)! W» 

we can obtain equation of motion for the GF as discussed 
in Sec. 2.1 (see Eqs. (8)— (11)). For the zero-order GF in 
MFA we get the following result: 

£(M) (kjW) = ^ + (g(k)-^)f 3 + A(k)f 1 
(e(k) - I A(k) | 2 

where the the holon dispersion and the gap function are 
given by the equations: 

e(k) = t(k) - 1 J2 - q) - J'(0)} K (55) 



(56) 



In comparison with MFA for RVB state considered pre- 
viously, in the GF approach we obtain additional renor- 
malization for the holon dispersion, Eq. (47), while the 
gap equation has the same form, Eq. (48), but with a 
different renormalized exchange interaction: 



= Jii{\-{S Z iS*)). 



(57) 



By using equation of motion method instead of MFA 
for RVB state we managed to take into account spinon 
correlations in the effective holon interaction (57) and 
obtained a simple formula for it by using the identi- 
ties: Sf S? = 0, Sr S+ = (1/2) - Sf in the last equa- 
tion. From Eq. (57) we get the following estimation for 
the effective exchange energy for the nearest neighbors: 
J' ~ (0.5 -7- 0.25) J if we assume the AFM Neel phase: 
(Sf Sj) = —1/4 or completely neglect AFM correlations: 
(Sf S?) = , respectively. 



A. Spectral density 

Spectral density for physical electron excitations in the 
lower Hubbard subband within the t-J model which is 
defined by the equation 



4,(k,w) - —Tm((XF\X?>)) u+ie , 
satisfies the sum rule: 

£ / du;A,(k, W ) = £<{X°",A7 }) 
= 2Q = 1 + S. 



(58) 



(59) 



In the MFA for the GF in terms of the HO, Eq. (21), we 
get: 

4(k,w)=Q Mk %-£ k ) +Qvk6(w + Eu), (60) 
where 

l n , (gk -fth l n (gk - A) \ 
k= 2^ — k= 2^ H£ — ' 

In MFA the spectral density (60) satisfies the sum 
rule (59) but have no incoherent background which ap- 
pears if one takes into account the self-energy corrections 
as shown in Ref. 36 . 

In the spinon-holon representation (29) the spectral 
density, as follows from Eq. (34), does not obey the sum 
rule (59): 

+ OC 

£ f duM^)(k,u,) = 

+ OC 

= £ J do; {-hm((ht b ™ \ h i } 



52{{hfb i<r ,h i bf}) = l. 



(62) 



MFA for the spinon-holon GF in Eq. (62) results in the 
spin-charge separation which defines the spectral den- 
sity (62) as a convolution of the anticommutator holon 
and the commutator spinon GF 28 : 



+oc />+oo 



„ J —OG J — OO 



dwida>2 

u> — u>i — U>2 



x [tanh(wi/2T) + coth(w 2 /2T)] 

x Im«A+ | h^U Im((,S£_ q | S£_ q )) W2 . (63) 

In this representation the spectral density of physical 
electron excitations is given by a superposition of the 
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spectral density for spinless fermion (holon) excitations 
and a background produced by spin excitations 28 . In 
MFA the holon spectral density in the paired state below 
T c is defined by the GF (54): 

= u k 5(uj + E(k)) + i\S(cu - 25(k)), (64) 

where E(k) = [(e(k) - ^ h ) 2 + | A(k) | 2 ] 1/2 is the quasi- 
particle energy and «k, £k are given by Eq. (61) but with 
quasiparticle energies for holons. The spectral density 
of electron excitations is measured in the angle-resolved 
photoemission (ARPES) experiments which provide in- 
formation concerning the symmetry of the gap function. 
As discussed below, the symmetry of superconducting 
gap function Ak defined by Eq. (27) appears to be d- 
wave, while the holon gap function A(k) in Eq. (51) 
has odd symmetry of the p-type as in the triplet pairing. 
This p-wave gap symmetry has never observed in ARPES 
experiments in cuprates. 

B. Gap symmetry 

For models with strong electron correlations the s-wave 
component of superconducting gap must be strongly sup- 
pressed due to on-site Coulomb correlations. For the t-J 
model it follows from the constraint of no double occu- 
pancy on a single site given by the identity: 

(CuaCi-a) = y^(Ck, CT C_ k ,- CT ) = 0, (65) 
k 

for the physical electron operators Ci i(J = Cj i(7 (l — ni_- a ). 
Since the anomalous correlation function (ck, CT c_k,-cr) is 
proportional to the gap function A(k), Eq. (65) imposes 
a certain constraint on the symmetry of the gap func- 
tion. In particular, for the solution (26) in MFA Eq. (65) 
reads 34 ' 35 : 

q 




For a tetragonal lattice the Fermi surface (FS) is invariant 
under the C4-axis rotation in the k-space and therefore 
to satisfy the condition (66) for E^ > the gap function 
Aq should change its sign along the FS. It means that 
the symmetric s-wave solution, A s (k x ,k y ) oc (cosq x + 
cosq y ) , does not fit Eq. (66), while the <i-wave solution 
with B lg or B 2g symmetry, Ad(k x ,k y ) = —Ad(k y ,k x ), 
satisfy the condition (66). In general, Eq. (66) should 
be considered as a constraint on the symmetry of the 
gap function in the superconducting phase and solutions 
which violate this constraint should be disregarded. 



Now we consider the symmetry of the spinon and holon 
order parameters. The condition (65) for the RVB order 
parameter reads 

B * = jf E w = a? E ^ i) F = °> ( 6? ) 

k k,q 

which also imposes a constraint on the symmetry of the 
spinon and holon order parameters. Since the symme- 
try of the holon pairing order parameter, Eq. (44), for 
spinless fermions is odd, imposed by their anticommuta- 
tion relations, it results in the odd symmetry of the gap 
A(k) in Eq. (48). For a tetragonal lattice the symme- 
try is given by the odd two-dimensional irreducible rep- 
resentation E u which can be modelled by the function: 
A(k) oc r]p (k) = (sink x ±smk y ) as in the p- wave triplet 
pairing. The same holds for the spinon order parameter 
inEq. (45), ^(k)cx 77±(k). 

The RVB order parameter B^ = (bij) = <PijF^ as a 
product of two antisymmetric order parameters with E u 
symmetry has either A (s-wave) or B\ g (d-wave) sym- 
metry. Namely, if we adopt only the nearest neighbor 
pairing for both the order parameters: 

<fij = <pi,i+a x {{Sj,i+a x — $j,i-a x ) ± {&j,i+a y — <5j',i-a y )}j 
^ij = Ftj,+a x {(fij,i+ax ~ fij,i—a x ) i {$j,i+a y ~ $j,i-a y )}, 

then for the Fourier transform of the RVB order param- 
eter we get 

s(k) = ^(k - q) ^ + (q) = 2b E < sin (^ - «*) 

q q 

± sin(fc y - q y )} {sinq x ± sing,,} 

= —b(cosk x ±cosk y ) (68) 

where b = (fii,i+a x F^ i+a ^ and the sign + (— ) in the 
last line corresponds to the same (different) signs in the 
first line. The solution (68) satisfies the condition (67) 
for both the s- and c?-wave symmetry: 5Z k B(k) oc 
^ k (cos k x ± cos k y ) = 0. 

However, the explicit solution for the holon order pa- 
rameter (50) shows its more complicated k-dependence 
which proves that the frequently used model for the near- 
est neighbor pairing is inadequate. For the solution (50) 
the condition (67) for the RVB order parameter reads: 

5>(k-q)F+(q) = 

k,q 

= Wk-q)^tanh^M = 0, (69) 
^ rK ' 2E(q) 2T 1 

k,q 

which also imposes a requirement on the symmetry of 
the singlet (RVB) order parameter which should have 
only the d-w&ve symmetry, B(k x , k y ) = —B(k y , k x ) , to 
satisfy this condition. 

It is interesting to compare results for the gap equa- 
tion derived in MFA for the t-J model and that one for 
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the Hubbard model. By applying the projection tech- 
nique for the GF for the Hubbard model one can get the 
following equation for the gap function (see, e.g., 39 ): 



A ija cx {XfNj) = (ci^Nj), 



(70) 



where we have used the identity for the Hubbard oper- 

From Eq. (70) follows 



ators, X- 



- V°J- Y^ 2 - n r- 



that the pairing occurs on one lattice site but in the dif- 
ferent Hubbard subbands. By using equation of motion 
for the GF ((X® 2 (t) \ Nj(t'))} the anomalous correlation 
function (Xf 2 Nj) can be calculated without any decou- 
pling that results in the same gap equation as in the MFA 
for the t-J model 39 : A ija = J tJ (X° a Xf)/Q where we 
have used the notation of the present paper. Therefore 
the symmetry constraint considered above is also applica- 
ble for the gap solutions obtained in the Hubbard model. 



C. Superconducting T c 

In the present section we compare the pairing tem- 
perature T c defined by the gap equation for physical 
electrons (27) and that one for holons, Eq. (51). We 
study T c as a function of corresponding chemical poten- 
tial v = [ijW for a given value of the pairing interac- 
tion A = J/W where all energies are measured in units 
of the renormalized half bandwidth for pairing fcrmions 
W = 4t: 



e(k) = #= 7 (k)+T 7 '(k), d(k)= A(k) 



W 



w 



(71) 



where r = t' eff /t eff . The gap Eqs. (27), (51) (or (56)) 
can be written in the form 



d(k) 



4A ^ 



7(k - q) d(q) 



2[(e(q)-^) 2 +d(q)]V2 



x tanh 



[(e(q)- I ,) 2 + rf(q)] 1 /2 



2T 



(72) 



Let us consider T c {v) for different pairing symmetry, 
d(k) = d a ?7„(k): 



s — wave : r/ s (k) = (cos k x + cos k y ), 
d — wave : ??d(k) = (cos k x — cos k v ) 1 
p — wave : r/^ (k) = (sin k x ± sin k y ) . 



(73) 



By integrating Eq. (72) with the corresponding symmetry 
parameter r] a (k) we obtain Eq. (72) for T = T c in the 
same form for any symmetry: 



1 



A N 
1 

_ 2 



2(e(k) 



tanh 



;(k) 



IT, 



(a) 



-N a (e) tanh- 



e — v 



IT, 



(a) ' 



(74) 



if we introduce an effective density of state (DOS) for the 
corresponding symmetry, a — s, d,p: 




FIG. 1. Effective density of states, Eq. (75), for s (solid 
line) , d (dot-dashed line) and p (dashed line) symmetry. 



k 



(75) 



which is normalized 

r+l i 

/ deN a (e) = 1, since — ^ (r/ a (k)) 2 = 1. 

k 

The results of calculation of the effective DOS, Eq. (75), 
for different symmetry a = s,d,p is presented on Fig. 1 
for t' = 0. From this dependence it is easy to draw a 
conclusion that the T c (v) function will follow the depen- 
dence N a (v) for corresponding symmetry since the effec- 
tive coupling constant V ~ XN a (v) reaches its maximum 
value at the maximum value of N a {v). In the logarith- 
mic approximation a solution for T c (p), Eq. (74), can be 
written in the conventional BCS form: 



T c ~ y/n(W-n)exp(-l/V), V-JN, 



Hi 



(76) 



but with the prefactor proportional to the Fermi energy, 
/i = Ef which can result in high-T c . The highest T c (u) 
appears for the cZ-wave pairing and the lowest T c (y) - for 
the p- wave pairing since in the former case the van Hove 
singularity gives a strong contribution for DOS ( [%(k = 
(±7r,0))] 2 = 4), while in the latter case it is completely 
suppressed ( [r/ p (k = (±7r, 0))] = ). This estimation can 
be checked by a direct numerical solution of Eq. (74) (see, 
e.g. 35 ' 36 and 49 for the p-wave symmetry). 



V. CONCLUSION 

In the present paper we consider superconducting pair- 
ing mediated by the exchange interaction which is generic 
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for system with strong electron correlations as cuprates. 
The mechanism of the exchange pairing is the lowering 
of kinetic energy of electron pairs due to their coherent 
hopping between different Hubbard subbands. Since the 
excitation energy of this hopping is much larger then the 
Fermi energy the retardation effects in the exchange in- 
teraction are negligible 39 which results in pairing of all 
electrons (holes) in the conduction band and a high-T c 
proportional to the Fermi energy, Eq. (76). 

To obtain estimation for the superconducting T c it 
is tempting to use a mean-field approximation for the 
exchange interaction within the t-J model. However, 
meaningful physical results can be obtained only if one 
takes into account strong electron correlations on a rig- 
orous basis which is provided by the Hubbard operator 
technique. Any auxiliary field representations applied 
within MFA inevitably violates rigorous commutation 
relations for HO which may result in unphysical con- 
clusions. In the present paper we have proved this by 
comparing the results for superconducting gap equation 
derived within the slave fermion - hard-core boson rep- 
resentation, Eq. (29), and HO technique for the Green 
functions (Sec. 2.1). In the former method the projected 
character of physical electron operators is neglected that 
results in double counting of empty states and violation of 
the sum rule, Eq. (62) . The spin-charge separation which 
occurs in MFA, Eqs. (40), (41), results in separate equa- 
tions for two order parameters for spinons and holons in- 
stead of one equation for physical electrons as in the GF 
method. The gap equation (51) for spinless fermions has 
only antisymmetric solutions which results in the p-wave 
gap for the quasiparticle excitations in superconducting 
state never observed in ARPES experiments in cuprates 
and much lower T c then for the <i-wave pairing given by 
Eq. (27) for physical electrons. 

The obtained results within MFA for the spinless 
fermion - hard core bosons in Sec. 2 appears to be identi- 
cal to the path-integral representation for the t-J model 
employed in Ref. 32 . Therefore the latter results have the 
same flaws as discussed above which casts doubts on ap- 
proaches based on the idea of spin-charge separation. Vi- 
olation of the local constraint in the two-band Hubbard 
model may lead also to unphysical results, as shown in 
Ref. 50 . 
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